Precise Measurement of the Absolute Yield of Fluorescence Photons in 
Atmospheric Gases 
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We have performed a measurement of the absolute yield of fluorescence photons at the Fermilab Test Beam. 
A systematic uncertainty at 5% level was achieved by the use of Cherenkov radiation as a reference calibration 
light source. A cross-check was performed by an independent calibration using a laser light source. A significant 
improvement on the energy scale uncertainty of Ultra-High Energy Cosmic Rays is expected. 



1. Introduction 

Fluorescence detection of Ultra High Energy 
Cosmic Rays (UHECRs) is a well established 
technique, pioneered by Fly's Eye pQ, and today 
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an integral part of the Pierre Auger i2; and the 
Telescope Array [3] experiments. Excitation of 
the atmospheric nitrogen by the charged parti- 
cles in the extensive air shower induces the emis- 
sion of fluorescence photons, mostly in the 300- 
400 nm range. A Fluorescence Detector (FD) 
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records this radiation to infer the cosmic ray en- 
ergy and the particle type. For this purpose, the 
fluorescence light yield from the charged parti- 
cles in the shower must be known for every emis- 
sion point along the shower axis. A correction 
is then applied to account for atmospheric effects 
between the shower and the telescope, enabling 
an accurate, quasi- calorimetric, primary energy 
determination. 

The uncertainty on the fluorescence light yield 
is one of the main systematic uncertainties on the 
cosmic ray energy determination by experiments 
that employ the fluorescence technique (e.g. 14% 
over a total 22% uncertainty for the Pierre Auger 
experiment). The AIRFLY collaboration has al- 
ready performed a very precise measurement of 
the fluorescence spectrum and its pressure de- 
pendence [4], as well as the dependence of the 
emission on the temperature and humidity [5]. 
AIRFLY measurements over electron kinetic en- 
ergies ranging from keV to GeV using several ac- 
celerators have also proven the proportionality of 
the fluorescence yield with the electron energy de- 
posit [6]. 

The final step in the precise characterization 
of the nitrogen fluorescence light emission is the 
measurement the absolute value of the yield for 
the main emission line at 337 nm. The AIRFLY 
strategy to reduce the systematic uncertainties is 
to calibrate the experimental apparatus in situ, 
using photons emitted by a well know process: 
Cherenkov radiation [7 . A second calibration 
method, with nearly independent systematic un- 
certainty, is based an absolutely calibrated laser 
light source. In this paper, we present prelimi- 
nary results from a series of dedicated measure- 
ments at the Fermi National Accelerator Labora- 
tory (Fermilab). 

2. Experimental method 

The measurements were performed at the Fer- 
milab Test Beam Facility. Most of the mea- 
surements were carried out using the 120 GeV 
proton beam, extracted from the Main Injector. 
Secondary beams of 32 GeV pions and 8 GeV 
positrons were also used. 

A sketch of the apparatus is shown in Fig. 




Veto counter 



Figure 1. Layout of the experimental apparatus 
used at the Fermilab Meson Test Beam 




Figure 2. AIRFLY chamber with the integrating 
sphere fitted inside. 
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Figure 3. The integrating sphere used in the mea- 
surements. The diffusive coating applied to the 
interior can be seen through two of the ports. 



[T] A fluorescence chamber made of a 3 mm 
thick stainless steel was placed in the beam path, 
with the corresponding flanges for windows, shut- 
ters, gauges, gas inlet and pump-out. Both the 
entrance and exit windows are 0.1 mm thick 
aluminum, with the entrance window before a 
18 cm long tube to provide additional length 
for Cherenkov light production. An integrating 
sphere was used to collect light produced inside 
the chamber (Fig. [5]). One of the ports of the 
sphere was open to a a gas-tight window fitted 
with a 337 nm filter and then coupled to the pho- 
ton detector. A Hamamatsu H7195P photomul- 
tiplier (PMT) tube, with good single photoelec- 
tron resolution, was used for photon detection. 
The optical field of view was defined by a 40 mm 
diameter acceptance cylinder placed between the 
integrating sphere's port and the filter, and by 
circular apertures of the same size placed in front 
of the PMT photocathode. A mechanical shutter 
remotely controlled allowed to take background 
measurements. The chamber was internally cov- 
ered with a black UV-absorbing material to avoid 
stray light. 

The purpose of the integrating sphere is 
twofold: it increases the light collection solid an- 



gle (and thus maximizes the signal to noise ratio) 
and it works as equalizer of the detection effi- 
ciency of the set-up for the isotropic fluorescence 
light and for the highly directional Cherenkov 
light. The integrating sphere (Fig. [3]) was built 
from two hollow aluminum half spheres coated 
with a material of very high diffusive reflectance. 
Light produced inside the sphere is collected over 
47r solid angle and at the same time isotropized by 
several diffusion bounces, so that the lambertian 
light output is independent of the original light 
distribution. A total of 4 ports were machined 
in the sphere, one was the detection port, two 
were for the beam entrance and exit, and the last 
one was opened at the top of the sphere. Two 
remotely controlled shutters could close the top 
and exit ports of the integrating sphere with plugs 
coated with the same diffusive material material 
using inside the sphere. With the exit port closed 
the Cherenkov is diffused back into the sphere and 
can reach the photon detector. With the exit port 
open the Cherenkov is absorbed by the chamber 
lining and thus only fluorescence can reach the 
photon detector. The top port compensates the 
open/closed position of the exit port to maintain 
always three open ports and thus the optical char- 
acteristics of the sphere, i 

The whole chamber was airtight and a remotely 
controlled system for gas and vacuum handling 
was used. Pure nitrogen and a dry-air mixture 
were used for the measurements and additionally 
helium and argon were used for background runs. 
The pressure, temperature and humidity inside 
the chamber were monitored remotely using the 
appropriate sensors. 

A set of particle counters was used for beam 
monitoring. At the entrance of the chamber a 10 
mm diameter finger counter was used for beam 
tagging. At the exit the beam was tagged with 
a Cherenkov counter, a 10 mm diameter cylindri- 
cal rod made of UV-transparent acrylic material. 
The rod was 30 mm long and allowed very good 
single particle resolution with fast timing. Before 
and after the chamber two big scintillator pads 
with a central 10 mm diameter hole for beam 
passage provided a veto for off-axis particles. 
The chamber and the counters were mounted on 
an optical breadboard for precise mounting and 
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alignment, and placed on a remotely movable ta- 
ble that allowed the alignment of the appara- 
tus with the beam by maximizing the finger and 
Cherenkov counters rates. The beam profile was 
monitored by wire chambers placed before and 
after the AIRFLY apparatus and was typically 3 
mm x 4 mm wide. 




Figure 4. The signal of a train of bunches for 
the Cherenkov rod (top plot) and the two veto 
counters (middle and bottom plots). Two par- 
ticles in two different bunches can be observed: 
the first one passes through the rod while the 
second one passes through both veto counters 
(and leaves also a small signal in the rod when 
it hits the PMT glass producing a small amount 
of Cherenkov in it). 



The trigger and DAQ were designed consider- 
ing the characteristics of the beam timing: within 



a 4 s beam spill, particles were grouped in trains 
of bunches. Trains were separated in time by 10 
/is, and bunches within a train were separated by 
19 ns. Typical conditions for data taking with 
the proton beam were 30 bunches per train, and 
a multiplicity of 2 ■ 10 5 particles per spill. The 
trigger logic was built from the coincidence of a 
train trigger gate, issued in correspondence of the 
arrival of each train of bunches, and a single par- 
ticle trigger gate, coming from a beam monitoring 
scintillator counter. Both triggers were provided 
by the Test Beam Facility. Whenever a trigger- 
was issued the signals from the scintillator coun- 
ters and photomultiplicrs were digitized by a 12- 
bit 500 MHz FADC and 600 samples (equivalent 
to 1.2 /us, containing the entire train of bunches) 
were saved in the FADC memory The data for 
the whole spill was stored in the FADC internal 
memory and then readout and saved to disk in 
about 40 s between spills. As an example, the 
ADC trace of the beam counters for one trigger 
is shown in Fig. [4] 

Data during the test beam was taken in a dif- 
ferent configurations. For each configuration runs 
of up to half hour (i.e. 30 spills) of data were 
acquired. Runs in the same conditions were re- 
peated periodically to improve statistics and to 
assure redundancy and consistency. 

3. Data analysis and results 

Offline data analysis for one run proceeds se- 
lecting single particles that cross the fiducial vol- 
ume of the camera, hence are tagged by the beam 
counters placed at the entrance and the exit of the 
chamber (Fig. [5]). For these selected particles we 
require that no signal is present in any of the veto 
counters. Additional analysis cuts are placed to 
discard trains that have unusually large number 
of particles passing through the veto counters, as 
this has been shown to improve the background 
conditions without excessively penalizing statis- 
tics. Once the clean single protons have been se- 
lected the PMT signal is analyzed and the pho- 
tons in coincidence with the selected protons are 
counted (Fig [6]). The signal S, in units of pho- 
toelectrons per beam particle, pbp, is then calcu- 
lated. 
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Figure 5. Spectrum of the Cherenkov rod. 
From left to right the three peaks correspond to 
pedestal, particles hitting the PMT window and 
single particles passing through the Cherenkov ra- 
diator. The fourth rightmost smaller peak corre- 
sponds to two particles in the same bunch. 



The measured signal in the a fluorescence mea- 
surement taken in a given gas, S 9 p^(meas), is 
given by: 

Sp^irneas) = Sp™ + Bp™, (1) 

where S 9 pl is the signal from the 337 nm band, 
and is attributable to background. The best 
determination of the overall background can be 
obtained combining measurements from different 
gases. From Eq. [IJ the difference of the measured 
signal in N2 and air is given by: 

AS FL = S*i - Sf L + B^ L - Bf L . (2) 

Since the beam related background and the pri- 
mary interactions are the same and secondary 
particle production is very similar in nitrogen and 
air backgrounds cancel in Eq. [2j Thus, 

AS Fi = ^l--L), (3) 

where tn 2 is the ratio of the 337 nm fluorescence 
in pure nitrogen to the signal in air. This ratio 
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Figure 6. PMT spectrum for a run. The events 
inside the single photoelectron peak are counted 
to define the signal. 



was measured by AIRFLY as a function of pres- 
sure in 0] and was cross-checked in our Fermilab 
experimental apparatus using a 241 Am radioac- 
tive source, giving 7.45 ± 0.10 at 1000 hPa. 

We measure AS FL = (16.83 ± 0.13) • 10~ 4 pbp. 
From Eq. G2 we derive the background-subtracted 
fluorescence signal: 

S^l = (19.44 ± 0.15) ■ 10" 4 pbp, (4) 

with the background (Eq. [lj estimated as a 3% 
of the signal. 

In the Cherenkov calibration measurement, 
both Cherenkov and fluorescence emission in the 
gas contribute to the measured signal: 

S c a ^(meas) = Sp™ + S c ™ + Bp™ + B C h, (5) 

where S c a H is the signal from Cherenkov light in 
the 337 nm band emitted in the gas under study, 
and Bch takes into account background origi- 
nating from the interaction of the beam particles 
in the exit port plug (that is independent of the 
gas filling the chamber). Bch can be estimated 
directly from the difference in the vacuum mea- 
surement in fluorescence and Cherenkov modes 

B C h = (2.63 ± 0.23) • 10~ 4 pbp, (6) 
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a value that is ~ 10% of the Cherenkov signal. 
We made several cross-check of this background, 
for example changing the exit port plug material 
to a thin Mylar foil where no light production is 
expected. 

We measure 5^(meas) = (32.89 ± 0.15) • 
10~ 4 pbp. Using Eq. [5]we obtain 

Sch = (10.27 ± 0.23) • 10~ 4 pbp, (7) 

and from it and Eq. |4] we derive 

c-N 2 

R N2 = = 1.893 ± 0.045, (8) 
"Cfl 

the ratio of fluorescence to Cherenkov 337 nm 
photons produced inside the chamber. 

In order to derive the absolute 337 nm yield 
in air, we performed a full Monte Carlo simula- 
tion of the setup, where all individual components 
are simulated according to measurements done in 
the laboratory. The absolute fluorescence yield in 
the simulation, Y^" c , which determines a corre- 
sponding expected Fluorescence/ Cherenkov ratio 
R^f c , can be scaled to match the measured value 
of Eq. [5] and obtain our measurement of the ab- 
solute fluorescence yield in air: 

Y air = 5.60 ± 0.13 7 337„m/MeV, (9) 

where the uncertainty is statical only. Our pre- 
liminary estimate for the systematic uncertainty 
of the measurement is better than 5%, dominated 
by the uncertainty in the wavelength dependence 
of the PMT quantum efficiency and the 337 filter 
transmission, accounting for the different spectral 
distributions of the fluorescence and Cherenkov 
signals. The detailed analysis of the uncertain- 
ties is still, however, ongoing. 

In the laser calibration method, we use a 
337 nm nitrogen laser and an NIST calibrated 
probe to determine the overall efficiency of the 
set-up, the number of detected photoelectrons per 
photon entering the sphere. Using this efficiency, 
the number in Eq. @] and a detailed Monte Carlo 
simulation of the laser calibration the value 

^LASER = 5.56 ± 0.0 7 7 337nm/MeV, (10) 

is obtained, which has an uncertainty nearly in- 
dependent of the one from the Cherenkov calibra- 
tion. The systematic uncertainty is dominated by 



the 5% uncertainty in the absolute calibration of 
the laser probe. 

4. Outlook 

The preliminary results on the absolute flu- 
orescence yield reported in this work are com- 
patible with the current values used in Fluores- 
cence Detector analysis (see Ref. [8] for a review 
of the different measurements.) We expect a fi- 
nal measurement with a systematic uncertainty 
below 5%, a significant improvement over pre- 
vious measurements, which will correspondingly 
improve the uncertainty on the energy scale of 
UHECR measurements. 

REFERENCES 

1. R. M. Baltrusaitis et al., Nucl. Instrum. Meth. 
A 240 (1985) 410. 

2. J. A. Abraham et al. [The Pierre Auger 
Collaboration], Nucl. Instrum. Meth. A 620 
(2010) 227 [arXiv:0907.4282l [astro-ph.IM]]. 

3. Y. Tameda et al, Nucl. Instrum. Meth. A 609 
(2009) 227. 

4. M. Ave et al. [AIRFLY Collabora- 
tion], Astropart. Phy s. 28 (2007) 41 
|arXiv:astro-ph/0703132| . 

5. M. Ave et al. [AIRFLY Collaboration], 
Nucl. Instrum. Meth. A 597 (2008) 50 
|arXiv:071 1.45831 [astro-ph]]. 

6. M. Ave et al. [AIRFLY Collaboration], Nucl. 
Instrum. Meth. A 597 (2008) 46. 

7. M. Bohacova et al. [AIRFLY Collabora- 
tion], Nucl. Instrum. Meth. A 597 (2008) 55 
|arXiv:0812.3649l [astro-ph]] . 

8. J. Rosado, F. Blanco and F. Arqueros, As- 
tropart. Phys. 34 (2010) 164 [arXiv: 1004.39711 
[astro-ph.IM]]. 



